Squeezed Light Generation in Nonlinear System with Chaotic Dynamics 
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The squeezing in a nonlinear system with chaotic dynamics is considered. The model describing 
interaction of collection of two-level àtoms with a single-mode self-consistent fïeld and an external 
field is analyzed. It is shown that in semiclassical limit, in contrast to the regular behaviour, the 
chaotic dynamics may result in: (i) an increase in squeezing, (ii) unstable squeezing and contraction 
of time intervals of squeezing on large enough times. The possibility of the experimental observation 
of the described effects is discussed. The obtained results are rather general and do not depend on 
the model under consideration. 

I. INTRODUCTION 



In recent years a great deal of interest has been focused on the investigation of chaotic dynamics in different 
physical processes [1,2]. Alongside with the study of dynamical chaos in the classical systems, in the middle of 70s 
morè and more researchers began to realize the importance of the problem of quantum chaos. The subject of quantum 
chaos is the study of some peculiarities in the behaviour of quantum systems which are classically chaotic [3-5]. A 
^ \ lot of interesting theoretical results have been obtained in the investigation of quantum chaos. One of the most 
important problem now is to study concrete physical systems with quantum chaos. Such systems have to meet certain 
requirements, i.e. they must be: a) semiclassical, b) chaotic in the classical limit, c) Hamiltonian (nondissipative). 
The last requirement is, apparently, not necessary but the study of the dissipative quantum chaos has only begun [6]. 
A very promising models in the study of quantum chaos are the models of quantum òptics. First of all, these are 
00 • models describing interaction between light and a collection of àtoms. Really, the semiclassical conditions and those of 
(semi) classical chaos can be easily met for such models (e.g., see review [7]). But the majority of the optical systems 
r"| ■ with classical chaos are dissipative (làsers, bistable devices, etc). Up to now the only exception was the BZT model 
i] [7,8] - the semiclassical Jaynes-Cummings model [9] without the assumption of the rotating-wave approximation 
+1» i (RWA). But the global chaos arises in the BZT model if only the RWA is violated [] ,and this requires very large 
atom densities and vàlues of dipolc moments in a experiments (numerical estimates see in [7,8]). Only recently 
there appeared some modifications of Jaynes-Cummings model assuming the Hamiltonian chaos in the semiclassical 
q-( description and within the framework of the RWA [10-13]. 

On the other hand, at present there are a number of purely quantum (nonclassical) effects arising in the interaction 
. £^ of light with àtoms (squeezing, photon antibunching, etc.) [14-17]. As a rule, these effects are more pronounced in the 
semiclassical limit. However, all the nonclassical effects with light have been studied for the integrable systems with 
regular behaviour 0. Therefore, in order to investigate the problems of quantum chaos, as well as those of quantum 
òptics, it is of interest to analyze the influence of the chaotic behaviour on the nonclassical optical effects. 

In this paper, a squeezed light generation in a nonlinear system with chaotic behaviour is considered. The model 
[12] is used as a concrete example, and the method of Í/N expansion into modification suggested in [18] is applied. 

The main result of the paper is the demonstration of the fact that in the semiclassical limit chaos may increase 
squeezing of light. A phenomenon of the unstable squeezing is described. The unstable character of squeezing is 
manifested in a great change of time intervals, during which the the squeezing is possible when the initial conditions 
or parameters of the systems change slightly. 

The paper isorganized as follows. Section 2 presents a brief description of the model [12] and the procedure of í/N 
expansion. Here also the equations are obtained which describe light squeezing in the semiclassical limit. Section 3 
deals with the conditions for transition to chaos in the model discussed. Section 4 is devoted to the comparison of 



1 The RWA is violated when the coupling constant between the àtoms and the field (in units of frequeney) becomes of the 
order of transition frequeney. 

2 The problem of a squeezing in chaotic system was mentioned in concluding remark of the paper [19]. But no analysis of the 
problem was given in [19]. 
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squeezing for regular and chaotic dynamics. In Discussion other models of quantum òptics are considered permitting 
semiclassical chaos and as consequence a possible enhance in squeezing. The conditions for observation of chaotic 
squeezing in experiments are also discussed. 

Prcliminary rcsults of the paper were presented in [20] . 

II. Í/N EXPANSION AND SQUEEZING 

Let a sample with a gas of N two-level systems (TLS) be placed in a ring single-mode, high-Q cavity of resonant 
frequency ui. Now let us consider the interaction between TV two-level àtoms, a self-consistent field and an external, 
classical, amplitude-modulatcd field of the form 

E cxt {t) = E F{t) cos ut (1) 

injected into the cavity in the plane z = 0. In eq. (1) Eq is the external field amplitude, and F(t) is a function, slowly 
changing in comparison with a carrier frequency u and dcfining the amplitude modulation. Using the approximation, 
known in the theory of optical bistability and làsers with an injected external signal as the "mean field model" [21], 
the Hamiltonian in the interaction representation can be written in the form 

H = NH N , H N = Hjc + fíext, (2) 
Hjc = iTiuj c (a + J- — aJ + ), (3) 
ff ext = ihu c GF(t)(a-a+), (4) 

[a,a+] = í/N, [J+, J_] = 2J Z /N, [J Z ,J ± ] = ±J ± /N, (5) 

where a + , a are the normalized operators of crcation and annihilation of photons in the cavity; J a = 
^ 1 EÍIi a a,! ( a = +i ~ > z ) are collective atòmic operators (a a .i are Pauli matrices describing the i-th atom); 
q = N/V is the dcnsity of TLS, d and u>o are the dipole matrix element and transition frequency, respectively; 
G = cuir/Iui 2 , l is characteristic sample size, c is light velocity , ur = dE /h is the Rabi frequency. In (2)-(4) the 
case of exact resonance between àtoms and field is considered ui — u)q. The atom- field coupling constant in frequency 
units lü c = (2irgd 2 u>o/h' 2 ) 1 /' 2 is called cooperative frequency [22] (in modern literature, u> c is also called "collective 
vacuum-field Rabi frequency" [23]). It characterizes the energy exchange from the àtoms to the field and back again. 
Such kind of oscillations have been observed experimentally in both optical and microwave domains [23-25]. 

One can easy see that eq. (3) represents the Hamiltonian of the Jaynes-Cammings model [9] , and eq. (4) belongs 
to the type "external source - classical current" , according to the classification suggested in [26] . 

Different modifications of driven Jaynes-Cammings model were introduced in the framework of investigation of the 
quantum field fluctuations including squeezing [27]. But this activity was limited only by the integrable case with 
regular dynamics. In contrast, the model (2)-(4) demonstrates both regular and chaotic dynamics. 

The quantum dynamics of our system is completely described by equation of motion for the density matrix 



ihp=[H,p] (6) 



Now we introducc new operators 



li 
X 1 = Rea = -(a + + a), X 2 = Iraa = -(a+ - a) (7) 

From the commutation rclationships (5), we obtain the uncertainty relation 

((5 Xl f)((8X 2 f) > J^, SX j =X j -(X j ) (8) 

The definition of squeezing has the form [16] 

S = 4N((SX 1 ) 2 ) < 1 (9) 

In order to study the dynamics of squeezing in the model (2)-(4) at N » 1, it is natural to use the method of Í/N 
expansion. It has been recently mentioned in [28] that the method can be applied not only to the integrable but also 
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to nonintegrable systems. In accordance with this method [29], the natural states for constructing the classical limit 
(N — > oo ) are generalized coherent states. In our case, the state | tjj) is the product of the Glauber coherent state 
a) [14] and of the spin coherent state | fi) [30-31]. For the model (2)-(4), the Í/N expansion is constructed in a 
manner analogous to the expansion for the Jaynes-Cummings model [18]. Therefore, l'd like to dwell here only on 
some principal moments. 

We define a generalized P-function as 

p(t)= f d 2 ad 2 f3P(a,í3,t) | I, (10) 



/r+oo ç+oo 
d 2 a = / d(Rea) / d(Ima) 
J — oc J — oo 



Using the commutation relations between the Bose operators and projectors [14] and the spin operators and projectors 
[31] 

a + | a)(a | = U a + \ a)(a |, 
a | a)(a \ =U a \a)(a\, 

\a)(a\a+ = U: + \a)(a\, (11) 
| a)(a | a = U* \ a)(a \, 



U a+ 



a 



u 

Nda' a 



a 



J+\[3){f3\ = U J+ \0)(0\, 
J-\(3)(Í3\ = Uj_ \f3){[3\, 
\f3){P\J + ^u: J _\f3)(P\, 
\(3){í3\J- = U* J+ \í3)(f3\, 



(12) 



Uj + = 



(3* 



Ndf3' 



Uj_ = 



1+I/3I 2 

from eq.(6), we have the equation of motion for the function P 



1+1/3 



0^d_ 

N d{3 



dP 

~dt 



dad(3 



(W aP P) 



c.c. 



(13) 



V a = 1+ P p + GF(t), V p = -(a + a*(3 2 ), W aP = -(3 2 /N (14) 

It follows from (13) and (14) that the variables a and j3 can be considered to be real. From eq. (13) one can obtain 
the following equations of motion for averages [18] 

^(ip)=u c {Vp), (15) 

^(Sipóíp/) = uj c {V v Sipf) + u) c (6ipV vl ) + u} c {W vv ,), (16) 

where Síp = ip — (ip), and (p is one of the variables a or (3. However, the equations (15) and (16) are not closed. So, 
let us use the Taylor's expansion of the function V v 
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and the analogous equation for W vv ,. Due to the uncertainty relation of the type (8), it is clear that S(p(t = 0) ~ N~ 1 / 2 . 
Therefore, substituting eq.(17) into (15) and (16) one can obtain a closed system of equations for observed vàlues in 
any order over X/N. In a zero order, the equation of motion takes the form 

j t { V >)=ü Jc {V tp ) {v) +0{l/N) (18) 

These equations are equivalent to the equations obtained from a coupled Maxwell-Bloch system in [12]. After replacing 
the variables j3 = tan(z/2) [30], eq. (18) can be written as follows 

x + lü 2 sinx = -2Guj 2 F(t), (19) 

where 

a=p/2, p= — —x, (J + ) = -~sinx, (J z ) = -i cosa; (20) 

Lú c Z Z 

It is known already [1,2,12] that eq. (19) assumes both regular and chaotic behaviour. The conditions for the transition 
to chaos in (18), (19) will be discussed in the next Section. 

In the first order in 1/N, one can obtain from eq. (16) with an account of (17) 

í((Ap) 2 ) = 2uj c cosx(A p Ax), 
dt 

^((Ax) 2 ) = -2u Jc (ApAx), (21) 

4:(ApAx) = lo c cosx((Ax) 2 ) - uj c ((Ap) 2 ), 
dt 

where 

((Ap) 2 ) = A((Sa) 2 ) + 2/N, ({Ax) 2 ) = 4R 2 ((S{3) 2 ) + 2/N, (22) 

(ApAx) = AR(SaS(3), R=(l + /3 2 ) -1 

In the present paper, we shall focuse our attention at the conditions when initially field is in the coherent state and 
the àtoms are in the ground state (J z = —1/2) corresponding to x(0) = and p(0) = po. The condition for squeezing 
in variables (22) takes the form 

S = N((A P ) 2 ) < 3 (23) 

Before we start investigating the nonlinear dynamics of the systems (19) and (21), let us discuss the condition for 
validity of the semislassical approach to squeezing. 

The time-scale for validity of semiclassical approach and 1 /N expansion in cooperative optical system were studied 
systematically in [32]. For regular dynamics, this time-scale has a power dependence on N. In contrast, it is 
proportional to log N for a chaotic dynamics [32] . 

Here we only estimate a required number of TLS, when the semiclassical description (18), (19) and (21) can be used 
to describe the dynamics of squeezing. Introduce the "convergence radius" 

d(t) = [((8a) 2 ) + ((SP) 2 )} 1/2 ~ [((Ap) 2 ) + ((Ax) 2 )} 1/2 « 1, (24) 

where d(0) ~ N^ 1 / 2 . If d(t) <C 1, then 1/7V expansion and consequently equations (18), (19) and (21) are correct. 
The behaviour of d(t) differs considerably for regular and chaotic dynamics. Equations (21), in fact, coincide^J with 
the equations arising in the definition of the maximal Lyapunov exponent [1]. It is known [1,2] that for the regular 
dynamics d(t) increase according to a power law and for the chaotic motion it is exponential 



3 The difference is in the fact that in the definition of Lyapunov exponent only the linearization near x(t), p(t) is considered 
and not the behaviour of the linear fluctuations as in (21). However, in this case the difference is insignificant. 
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d(t) = d(0)exp(\Lü c t), 



(25) 



where A > is the maximal Lyapunov exponent. The dependence (25) is related to the presence of the strong local 
instability of the chaotic motion. We will consider the effective squeezing to occur for the time t\ , being of the order 
of several uj^ 1 , and Í/N cxpansion is correct, if d < 0.01 . With this assumption, from (25) it follows that one need 
to have N ~ 10 5 -j- 10 6 for correct semiclassical description of squeezing at chaos conditions. This simple estimation 
is in a good agreement with the results of numerical simulation [20]. 



III. CHAOS 



In this section, following [12], we will briefly discuss the conditions for transition to chaos in the model (18), (19). 
We will need this information below. 

The motion of a pendulum without perturbation (G = 0) is periòdic, and in phase plane it has two types of the 
fixed points: elliptic, with coordinates p — 0, x — 2nn (n = 0, ±1,±2, ...) corresponding (see (20)) to the initial 
populations of de-excited levels of TLS, and hyperbolic, with coordinates p — 0, x — 7r(2n + 1) (n = 0, ±1, ...) 
corresponding to the complete filling of the upper levels of TLS. The pendulum separatrix (a special trajectory in the 
phase plane, separating the vibrational and rotational motion and going through the hyperbolic points) corresponds 
to the complete energy transformation from the àtoms to the field and back again. 

The dynamics of (18), (19) at G ^ depends considerably on the number of harmonies in the spectrum of F(t). 
Let F(t) = sinílí. The criteria for the transition to chaos obtained by means of the Chirikov's resonance overlapping 
mcthod [1-4] are different in two limiting cases: k » 1 and k -C 1, k = 2Gui^/iï 2 . 

1) If k <C 1 and for fi > oj c , a narrow stochastic layer is found in the vicinity of the separatrix; the rest of the phase 
space is filled by periòdic trajectòries. At Q < lu c in the vicinity of the separatrix there appears a broad stochastic 
layer which occupies the major part of the phase space, except the area in the vicinity of the elliptic point. 

2) If k ^> 1, the criterion for the transition to chaos 

K > 1, K = const^yj, const ~ 10 (26) 
can be easily satisfied. The oscillation amplitude in this case may be rather large 

| Pmax N = —FT~ ( 27 ) 

According to eq. (26) at — > chaos is always present, but in this case the diffusion rate falls in proportional to íí, 
and the chaotic motion in the system is slow (adiabatic chaos) . Of great importance is the fact that in the case under 
consideration the transition to chaos is possible, if only the lower levels of TLS are initially occupied. 

A numerical analysis shows that at k ~ 1 the dynamics of the system is qualitatively analogous to the case of 

K > 1. 

When the number of harmonies for F(t) is increased, the chaotic properties of the system become strongly. 

Fig. 1 (a,b) and 2 show examples of chaotic, regular oscillations and adiabatic chaos, respectively (r = üj~H). Fig 
3 demonstrates the presence of the local instability in the chaotic motion (a) and its absence in the regular motion 
(b). 

Dissipative analog of such model was considered in [19] (both atòmic and field dumping were ineluded). In that 
case chaos is transient. 



IV. SQUEEZING AND CHAOS 

Let us consider now the dynamics of light squeezing. As the system (18), (19) is Hamiltonian, the Liouville theorem 
is always true. However, the change in the shapc of phase volume for regular and chaotic dynamics is different. It 
is known [1,2] that for chaos the change of the phase volume shape is the strongest and fastest. This is related to 
the fact that in the Hamiltonian systems the mechanisms for the onset of chaos are the phase contour stretching in 
one direction (and squeezing in another) and its further folding. This procedure of squeezing and folding is repeated 
several times. Eventually, the structure of the phase volume contour becomes complicated, and its perimeter length 
increases exponentially with time. The time intervals, when full squeezing is possible, become very short and alternate 
irregularly due to the chaotic motion. The dynamics of phase volume contour for chaos is very sensitive to any changes 
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in the initial conditions and the system's parameters. This peculiarity of the chaotic system is naturally manifested 
in the behaviour of the system (21) determining the dispersion behaviour. 

Though in the evolution process the phase volume shape of a nonlinear system with regular dynamics also undergoes 
considerable dcformation, the stretching is less and rathcr regular. The increase in length of the phase volume contour 
has a power dependence. Therefore, the time intervals of squeezing have to be long even for large times. The regular 
character of motion has also to result in a weak dependence of the intervals on the parameters and initial conditions. 

These simple qualitative considerations have been verified numcrically. The squeezing (23) has been determined by 
solving the system (19) and (21) when F(t) = sinfií. The Rungc-Kutta method of the fourth order has been applicd 
with an accuracy up to 10~ 9 -j- 10~ 13 . The precision of these calculations has been monitored by checking that the 
time-invariant of motion 

L = p/2 - cos x + 2xGsintp + ÏII, (28) 



dx dL dp dL dip ^ dl dL 

dr dp' dr dx ' dr dr dip' ' c 

is satisfied. Fig 4a shows the largest squeezing S^in for the time r = IOw" 1 as a function of the external perturbation 
G, and Fig 5a shows the same as a function of the external field frequency Í2. The increase of d during the same 
time characterizes a degree of the dynamical instability (Figs. 4b and 5b). The conclusion about the character of 
oscillations in the system (regular (R), chaotic (C), adiabatic chaos (AC)) have been made on the basis of a phase 
portrait and the behaviour of d(t) for times r = 200. A boundary in the space of parameters between chaos and 
adiabatic chaos is, naturally, quite relative. It can be easily seen from figures that the degree of squeezing is larger 
for chaos. The value of squeezing in this case makes up 10~ 2 -j- 10 ~ 3 . 

Several comparatively low vàlues of squeezing (S^in — 0.1) observed under chaos for some parameter vàlues have 
been due to the following rcasons: 

1) Weak statistical properties of the chaotic oscillations ( e.g., point 1 in figures 4 and 5) which are, as rule, related 
to the fact that these parameters are in the vicinity of the area of parameters corresponding to the regular dynamics. 

2) Adiabatic chaotic dynamics. A characteristic period of the oscillations in this case (see Fig. 2) is compared with 
the time during which the squeezing was determined (ti = 10). 

3) Anomalously strong instability of motion for r < 10 (see point 2 in figure 4a, b). The analysis of the behaviour 
S(t) at this parameter value shows that after a short fall S(t) quickly increases during the interval (ti = 10). This is 
due to fast and strong deformation of the phase volume. However, strong squeezing would be also possible further, 
e.g. for n = 20 S min ~ 2.6 • 10~ 3 . 

Fig. 6 shows temporary intervals of squeezing (S < 3), when the dynamics of the system is regular or chaotic. It is 
seen that the intervals are quite long and regular even at r > 10 for the regular dynamics (a) and short and irregular 
for the chaotic behaviour (d). It is noteworthy that at r > 30 the intervals of squeezing become very short and 
their finding in numerical calculation makes difficultics. The time interval of squeezing for the chaotic motion is more 
sensitive to slight changes of the initial conditions than for the regular one (compare (b) and (e), (c) and (f)). 

Thus, the chaotic dynamics leads to the short, irregular and unstable squeezing at r > 10. It should be noted that 
for N 3> 1 the value of the maximal squeezing does not, in fact depend on N. 



V. DISCUSSION 



Thus, the main result of the paper is somewhat paradoxical: a nonlinear system, which is in fact a generator of 
noise, may be more effective in suppressing quantum noise. The maximal squeezing for the chaotic motion may be 
by a factor of 10 3 more than for the regular dynamics. 

We believe that the effect discussed may be observed in a modification of the experiments on squeezing, when the 
light interacts with an ensemble of two-level àtoms [33]. The amplitude modulation of the external field leading to 
chaos may be not only sinusoidal but may also represent a sequence of short light pulses. Such a sequence could be 
generated by a mode locked làser. 

It should be also noted that in this paper a ring cavity is considered. But it can be shown that the transition to 
chaos is also possible for other cavity configurations. For example, in the case of a Fabry-Perot cavity eq. (19) has to 
be replaced by the pendulum type equation with a Bessel function instead of sinus (see the second reference in [12]). 
Such an equation has also chaotic solutions, but Hamiltonian chaos is transient [12]. 
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Bcsides our model, an analogous scheme can be applied to other models of cooperative optical system with chaotic 
dynamics in the semiclassical limit: 

1) The BZT model [8] and also its generalization for the multi-level systems [34]. In these models the global chaos is 
possible if only the RWA is violated. 

2) The Shepelyansky model [10]: an ensemble of 3-level systems interacting with two modes of self-consistent field. 
The global chaos is possible within the framework of the RWA when the dipole moments of one- and two photon 
transitions are commensurable. 

3) The ensemble of TLS interacting with the self-consistent field and the resonance external field of constant amplitude 
[11,13]. The transition to chaos is also possible within the framework of the RWA. According to classification of [26] 
this model belongs to the type "external source - classical field" . 

Consequently, in all the cases one should expect an increase in squeezing when the conditions for chaos are satisfied. 

Another class of models assuming the 1/N expansion and chaotic dynamics are those models which describe the 
parametric interaction of light waves. This can be easily understood if to take into account the fact that any cooperative 
òptic process may correspond to a parametric process (e.g., see [35]). At present time, the several situations are known 
when múltiple parametric interaction of both sealar light waves [36] and waves of different polarizations [37] appears 
at chaotic spatial evolution. 

The consideration of squeezing at chaotic light waves evolution will be the subject of our future publication. 
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Fig. 1. Nonlinear dynamics: (a) chaotic behaviour (G = 1., 
{}/(ü c = 1.); (b) regular behaviour (G = 0.1,/2/eu c = 1.). The 
initial conditions are i(0) = p(0) = 0. 
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Fig. 2. Adiabatic chaos: G = 0.5, .ft/wc = 0.2. The initial condi- 
tions are the same. 
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Fig. 3. (a) Local instability for chaotic behaviour shown in Fig. 
la; (b) the absence of local instability for regular motion. The 
parameter vàlues are the same as in Fig. lb. Tn both cases N = 10 6 . 
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Fig. 4. Maximal squeezing S min and the value of dynamical insta- 
bility In d for the time r = 10 as a function of the external per- 
turbation G. The initial conditions: x(0) = p(0) = 0.; N = 10 6 , 
/2/oí c = 1 . Point 1 corresponds to G = 0. 1 , and point 2 to G = 1.1. 
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Fig. 5. MaximalsqueezingS min (r = 10) (a) and the corresponding 
value of the dynamical instability \nd{r= 10) (b) as a function 
of í2/ü) c , The initial conditions are the same; G = 0.5, N = 10 10 . 
Point 1 corresponds to í)/w c = 1.1. 
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Fig. 6. Time intervals of squeezing (dark) for regular (a)-(c) and chaotic (d)-(f) dynamics. The parameter vàlues in figures (a)-(c) 
are the same as in Fig.lb, and in figures (d)-(f) they are the same as in Fig.la. For all the figures N = 10 10 , jr(0) = 0., p(0) = S. 
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